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Abstract 

In recent years, several fractional generalizations of the usual Kramers- 
Fokker-Planck equation have been presented. Using an idea of Fogedby 
[H.C. Fogedby, Phys. Rev. E 50, 041103 (1994), we show how these 
equations are related to Langevin equations via the procedure of subor- 
dination. 

Introduction. - Some 70 years ago, Kramers [1; considered the motion of a 
Brownian particle subject to a space-dependent force F(x) per unit mass. His 
goal was to compute the joint probability distribution /(x, u, t) for finding a 
particle at time t at the position x with the velocity u. For this quantity he 
could derive the famous Kramers- Fokker-Planck (KFP) equation [51 [3] 



d 

— + u • V x + F(x) • V, 
at 



/(x,u,f) =£ FP /(x,u,t) (1) 



where £fp is the Fokker-Planck collision operator 

£ F pf = lVu-(uf)+DA u f. (2) 

As is well known, Eq. |T]) corresponds to the Langevin equations (see, e.g., 
Ref. 0) 

^x(i)=u(t), ^u(t)=F(x)- 7 u(t) + r(t) (3) 

where T(t) obeys white noise statistics. It describes a Brownian particle which is 
subject to the relation (x 2 ) ~ Dt where D is the diffusion coefficient. In many 
complex systems, this relation is violated, however. In fact, one often finds 
(x 2 ) ~ D a t a with a/1 which is described as "anomalous" or "strange" diffu- 
sion. Here, D a is a generalized diffusion coefficient with units [D a ] = m 2 s~ a . 
Depending on a, such a process is called subdiffusive (a < 1), superdiffusive 
(1 < a < 2), ballistic (a = 2), or turbulent-diffusive (a = 3). 
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As was shown by several authors, strange diffusion may be described by 
fractional generalizations of Eq. (JT|) (for a review including discussions of various 
applications see e.g., Ref. [1]). However, the latter may differ in the way the 
fractional character is incorporated. Thus it comes as no surprise that three 
different types of fractional KFP equations may be found in the literature. E.g., 
Metzler and Klafter [51 [7] proposed the equation 



d/(x,u,t) 
dt 



[-u • V, - F(x) • V„ + C FP ] ls D\- 6 f{^, u, t) 



(4) 



which they obtained by means of a non-Markovian generalization of the Chapman- 
Kolmogorov equation. Another kind of fractional KFP equation has been pro- 
posed by Barkai and Silbey [5], namely 



d_ 

dt 



■u-V x + F(x) • V,, 



/(x,u,t) =£ FP75 A 1 " 5 /(x,u,i) 



(5) 



and D] 5 is the 



a-<5 



where 75 is a damping coefficient whose units are [75J = s" 
fractional time derivative whose Laplace space representation reads D\ i-r 
A 1-5 . Finally, employing the concept of continuous time random walks (CTRWs), 
Friedrich and co-workers [12 [TO] were able to derive the equation 



— + u • V x + F(x) • V, 



/(x,u,t) =£ FP75 2? t 1 - 5 /(x,u,i) 



(6) 



Here, T>\ s denotes a fractional substantial derivative which can be written in 



n-8 



Laplace space as V t ~ <-> [A + u • V x + F(x) • V„] 

In this letter we address the important question of how these three fractional 
KFP equations are connected to each other. As is well known, fractional dif- 
fusion equations can be linked to CTRWs. [TT] And according to Fogedby [12], 
the latter can in turn be linked to sets of Langevin equations. This fact will 
be exploited below in order to gain insight into the nature of the stochastic 
processes underlying the three scenarios that were just described. For the sake 
of simplicity, we restrict ourselves to one-dimensional problems with no external 
force. 

Langevin approach to fractional diffusion equations. - In the spirit of Fogedby, 
let us first consider a stochastic process which is described by the following sys- 
tem of Langevin equations: 



d 



u(s) = — ju(s) + r(s) . 



■t(s)=r 1 (s) 



(7) 



ds ds 
Here, the variable s is to be interpreted as an internal time, whereas t is the 
physical (wall-clock) time. Moreover, T and 77 are described, respectively, by 
Gaussian and one-sided Levy stable distributions [denoted by L$(x)} [13]. Math- 
ematically speaking, the stochastic process u(s) is subordinated by the t(s) pro- 
cess. The latter is invertible, and the probability density of finding the internal 
time s at time t is given by 

'' ' L s (t/M 1/S ) 



p(s,t) oc 



els 



1 



(8) 



which is called the inverse one-sided Levy stable distribution. It is a solution of 
the equation 

-p(s,t) = -— ls Dl- 5 p(s,t) , (9) 



2 



and its Laplace transform reads p(s, A) oc A 5-1 exp(— js^ S s). Assuming that the 
stochastic processes t(s) and u(s) are statistically independent, the probability 
P(u, t) of finding the velocity u at time t can be written as 



P(u,t) 



Po(u, s)p(s, t) ds 



(10) 



where the distribution function Pq(u, s) is a solution of the standard diffusion 
equation, 



dP (u,s) 
ds 



= C FP P (u,s) 



(11) 



From Eqs. l[9])-(|lip it then follows that P(u,t) satisfies the fractional diffusion 
equation 



^=W>J-'P(u,t). 



(12) 



The idea of representing the solution of a fractional diffusion equation like 
Eq. (|12p as a superposition of Gaussians goes back to Barkai [H] . In the follow- 
ing, we will extend this method from velocity space to phase (position-velocity) 
space. 

The fractional KFP equation by Metzler and Klafter. - Let us now consider 
the Langevin system 

■jj-x(s) =u(s), -^j-u(s) = --yu(s) + F(s) , ^-t(s)=r](s) (13) 

which is closely related to Eq. ©. Here, both x(s) and u(s) are subordinated 
by the same t(s) process. In analogy with Eq. (fTX)]) , the probability distribution 
f(x, u, t) can be written as 



f(x,u,t) 



fo(x,u,s)p(s,t) ds 



where fo(x,u, s) is the solution of the Fokker-Planck equation 

d 



d_ 

ds 



dx 



f (x, u, s) = C F pf (x, u, s) 



Using these relations together with Eq. |9]), one obtains 



df{x,u,t) 
dt 



d_ 

dx 



-FP 



1S D\- S f(x,u,t) . 



(14) 



(15) 



(16) 



This is the fractional generalization of the usual KFP equation considered by 
Metzler and Klafter [SJ El [Z] . Thus we have shown that the corresponding 
stochastic process is given by Eq. (TIB"]) . 

The fractional KFP equation by Barkai and Silbey. - Next, we want to 
consider the Langevin system 



dt 



x{t) = u{t) 



ds 



u(s) = — ju(s) + r(s) 



ds 



t(s) = T)(s) 



(17) 
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Here, the velocity coordinate is subordinated by the internal time process while 
the evolution of the space variable is in physical time. These equations may be 
rewritten as 



ds 



x(s) = u(s)r)(s) , 



^-u(s) = - 1 u(s)+r(s), ±t(s)=r)(s). (18) 
as as 



For any specific realization of ?y(s), one can view this as a stochastic process 
which only depends on the Gaussian variable T(s). The corresponding proba- 
bility distribution fo(x,u,s) is subject to the KFP-type equation 



8 . . d 



f (x, u, s) = C F pfo(x, u, s) . 



(19) 



The solution of this equation is a Gaussian probability distribution with the 
second order moments defined by 



as 

— {xu)(s) 
as 



-2 7 (u 2 )(s) + 2D, 

-7 (xu)(s) +r](s) (u 2 )(s), 



^{x 2 )(s) = 2 V (s) (xu)(s) . 
as 

For simplicity, we first consider the case 7 = in which one obtains 
(u 2 )(s(t)) = 2Ds(t), 
j t {xu)(s(t)) = {u 2 )(s(t)), 

j t {x 2 ){s{t)) = 2{xu)(s(t)). 

Introducing the auxiliary variables a(t) and £(£) via 



(20) 



(21) 



dt 



a(t) = s(t) . 



dt 



E(t) = a(t) , 



finds 



(u 2 )(s{t)) = 2Ds(t) , (xu)(s(t)) = 2Da(t), (x 2 ){s(t)) = 4D £(<) . 
and the characteristic function 



Z(k,a,-) — J dx J du f(x,u, •) exp[ifca; + iau] 

of fo is obtained as 

Z Q (k, a, s, a, £) = exp [-Da 2 s(t) - 2Dak a(t) - 2Dk 2 E(t)] 



(22) 
(23) 
(24) 
(25) 



[We note that the latter can also be calculated directly from Eq. (JT7J) .] 

Now, the stochastic process i](s) defines a probability distribution W(s, a, £, t) 
where s(t), <r(t), and £(i) are related to r/(s) via Eqs. (fT5|) and ([2"2"|) . We assume 
that this function satisfies the equation 



Odd 



W(s, a, E, t) = -— ls D]- s W{s, a, E, t) (26) 
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which is the natural generalization of Eq. ©. In analogy with Eq. (fl"T)|) , the 
generic characteristic function of f(x,u,t) can thus be written as 



Z(k,a,t) — J ds J da J <i£ Zo(k,a,s,a,'S)W(s,a,'E,t) . 



(27) 



It is straightforward to show that the corresponding distribution function f(x, u, t) 
obeys the fractional KFP equation 



8 



8 



dt dx 



f(x, u, t) = C FP j s D] s f(x, u, t) 



(28) 



with 7 = 0. 

The case 7 ^ is only slightly more difficult. Eq. (|28|) yields 



— -fc — 

dt da 



Z(k, a, t) 



Using Eq. ([25]) . we note that 



da 



and 



9 2 

7a- h Ua 

oa 



Z 



Zn = 



ja-^- + Da 2 
oa 



ls D]- 5 Z{k,a,t). 



Z n 



J d 
2 7 s - 1)— +70-— 
as da 



(29) 

(30) 
(31) 



By means of Eq. ([2"?|) . we then obtain 
d 



d d 
dl +(T d^ + S d-a 



W{s,a,T,,t) 



9 , n v d 
— 2 7 s-l) + 7— a 
ds da 



ls Dl- 5 W(s,a,Z,t) 
(32) 

which is a generalization of Eq. ©. So, starting with Eq. (JT7J) and assuming 
that the stochastic process r](s) defines a probability distribution W(s, a, S, t) 
which satisfies Eq. (|3"2"]) . the fractional KFP equation a la Barkai and Silbcy 
[Eq. ((2HJ)] holds. 

The fractional KFP equation by Friedrich et al. - Finally we address the 
generalized KFP equation with retardation proposed in Ref. [5]. It reads 







d 



dt dx 



f(x, u, t) = £ F p<y s T>} s f(x, u, t) 



(33) 



where T>\ s is the fractional substantial derivative introduced above. According 
to Ref. [TU], this equation can also be written as 



d 
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dt U dx 



f(x, u, t) = C FP / Qit-t^e^ 1 - 1 '^ 9 * f(x,u,t')dt' (34) 



if the memory kernel Q(t — t') is chosen appropiately. In this case, the charac- 
teristic function satisfies the equation 



— -k— 

dt da 



Z(k, a, t) = 



Q{t-t') 



ja^- + Da 2 
da 



Z(k,a + k(t -t'),t')d£ . 



(35) 
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and one finds 

Zo(k, a + k(t — t'), s, a, E) = exp 



-Da 2 s - 2Daka - 2Dk 2 Y, 



(36) 



where we have introduced the new variables a — a + s(t — t') and E = E + er(i — 
t') + s(t — t') 2 /2. Using relations similar to those in Eqs. (|3"0|) and (j3~lj) as well 
as the ansatz (|27l) . we obtain the evolution equation 



d d d 



Q{t-t') 



W(s,o-,E,t) 

— (27S- 1) +7TpCT 

as acr 



xW(s, a-s(t- t'), £ - cr(t - t') + -{t - t'f,t') dt' 



(37) 



for W(s, a, E, t). The latter is simply a retarded version of Eq. (I32|) . To clarify 
the difference between Eq. ff32]) and Eq. (l37l) we introduce the shifted variables 



si and E = E — ert + ft 2 in the sense that 



W(s, a, E, t) = FU(s, cr - si, E - at + -i 2 , i) 



(38) 



holds. Consequently Eq. ((37|) can be written as 



Q(t-0 



si) 



W/(s,cr,E,i')di'. 

(39) 

The difference in the two approaches of Barkai and Silbey and Friedrich et 
al. can be traced back to the Langevin process (|22p. Integration yields for 
instance 



a(t) 



dt's{t') 



<(*) 



ds's'r)(s') 



(40) 



where 77(5) is a stochastic process. Since one has to integrate over the product 
s?7(s) a stochastic interpretation of this integral is needed. A similar situa- 
tion arises in stochastic processes involving multiplicative white noise sources, 
where different interpretations of similar integrals have been given by Ito and 
Stratonovich (for a discussion see e.g. [2])- We conjecture that different inter- 
pretations of this integral and the integral arising for the variable E(t) lies at the 
origin of the two different fractional equations for the probability distributions 
W(s,cr,E,i). 

Conclusions. - Three different types of fractional generalizations of the 
KFP equation describing anomalous diffusion of inertial particles can be found 
in the literature. Based on the idea of subordination, which is equivalent to 
the introduction of an intrinsic, fluctuating time, we have clarified the meaning 
of these different equations. Whereas in the approach of Metzler and Klafter 
[5] - [7] both position and velocity depend on the intrinsic time, the approach 
of Barkai and Silbey [8] and Friedrich et. al. [9], [10] assumes that only the 
velocity is subjected to the subordination procedure. We conjecture that the 
difference between the approaches of Barkai and Silbey and Friedrich et al. is 
due to different interpretations of stochastic integrals. 
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